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The nature of cosmological solutions for a homogeneous, anisotropic Universe given by 
a Bianchi type-I (BI) model in the presence of a Cosmological constant A is investigated by 
Tj- ■ taking into account dissipative process due to viscosity. The system in question is thoroughly 

studied both analytically and numerically. It is shown the viscosity, as well as the A term 
! exhibit essential influence on the character of the solutions. In particular a negative A gives 

rise to an ever-expanding Universe, whereas, a suitable choice of initial conditions plus a 
positive A can result in a singularity-free oscillatory mode of expansion. For some special 
cases it is possible to obtain oscillations in the exponential mode of expansion of the BI 
model even with a negative A, where oscillations arise by virtue of viscosity. 



PACS numbers: 03.65.Pm, 04.20.Jb and 04.20.Ha 

Keywords: Bianchi type I (BI) model, Cosmological constant, viscous fluid 



> 

in 
o 
o 

O I. INTRODUCTION 

"o 

The investigation of relativistic cosmological models usually has the energy momentum tensor 
of matter generated by a perfect fluid. To consider more realistic models one must take into 
account the viscosity mechanisms, which have already attracted the attention of many researchers. 
Misner (J, 0] suggested that strong dissipative due to the neutrino viscosity may considerably 
^ ■ reduce the anisotropy of the blackbody radiation. Viscosity mechanism in cosmology can explain 
the anomalously high entropy per baryon in the present universe ||3|,|4|]. Bulk viscosity associated 
with the grand-unified-theory phase transition o may lead to an inflationary scenario |6j|2ll3]. 

A uniform cosmological model filled with fluid which possesses pressure and second (bulk) 
viscosity was developed by Murphy [9]. The solutions that he found exhibit an interesting feature 
that the big bang type singularity appears in the infinite past. Exact solutions of the isotropic 
homogeneous cosmology for open, closed and flat universe have been found by Santos et al [ 10], 
with the bulk viscosity being a power function of energy density. 

The nature of cosmological solutions for homogeneous Bianchi type I (BI) model was inves- 
tigated by Belinsky and Khalatnikov Eli by taking into account dissipative process due to vis- 
cosity. They showed that viscosity cannot remove the cosmological singularity but results in a 
qualitatively new behavior of the solutions near singularity. They found the remarkable property 
that during the time of the big bang matter is created by the gravitational field. BI solutions in case 
of stiff matter with a shear viscosity being the power function of of energy density were obtained 
by Banerjee II 1211 . whereas BI models with bulk viscosity (rj) that is a power function of energy 
density e and when the universe is filled with stiff matter were studied by Huang lfl3ll . The effect 
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of bulk viscosity, with a time varying bulk viscous coefficient, on the evolution of isotropic FRW 
models was investigated in the context of open thermodynamics system was studied by Desikan 
fLUl . This study was further developed by Krori and Mukherjee lfl5ll for anisotropic Bianchi mod- 
els. Cosmological solutions with nonlinear bulk viscosity were obtained in lfl6ll . Models with both 
shear and bulk viscosity were investigated in 001. 

Though Murphy [9] claimed that the introduction of bulk viscosity can avoid the initial singu- 
larity at finite past, results obtained in 11911 show that, it is, in general, not valid, since for some 
cases big bang singularity occurs in finite past. 

We studied a self-consistent system of the nonlinear spinor and/or scalar fields in a BI spacetime 
in presence of a perfect fluid and a A term lEolElll in order to clarify whether the presence of a sin- 
gular point an inherent property of the relativistic cosmological models or is it only a consequence 
of specific simplifying assumptions underlying these models? Recently we have considered a sys- 
tem of nonlinear spinor field in a BI Universe filled with viscous fluid Il22ll . Since the viscous fluid 
itself presents a growing interest, we have studied the influence of viscous fluid and A term in the 
evolution of the BI Universe Il23ll . In that paper we consider only some special cases those allow 
exact solutions. In this paper along with those special cases we study some general cases, giving a 
qualitative analysis of the system of equations. We also perform some numerical calculations and 
compare the results obtained with those given in some pioneering papers in this field, e.g. Ill ill 



H. DERIVATION OF BASIC EQUATIONS 

Using the variational principle in this section we derive the fundamental equations for the grav- 
itational field from the action (12.11) : 



y(g;e) = J J^V^dQ. (2.1) 

with 

=Sf = =2grav. ~\~ (2.2) 

The gravitational part of the Lagrangian (12.21) S? gmv . is given by a Bianchi type-I metric, 
whereas the term Jzf v f describes a viscous fluid. 

We also write the expressions for the metric functions explicitly in terms of the volume scale 
T defined bellow (I2.18I) . Defining Hubble constant (12.281) in analogy with a flat Friedmann- 
Robertson-Walker (FRW) Universe, we also derive the system of equations for T, H and £, with e 
being the energy density of the viscous fluid, which plays the central role here. 



A. The gravitational field 

As a gravitational field we consider the Bianchi type I (BI) cosmological model. It is the sim- 
plest model of anisotropic universe that describes a homogeneous and spatially flat space-time 
and if filled with perfect fluid with the equation of state p = £e, £ < 1, it eventually evolves 
into a FRW universe ll24ll . The isotropy of present-day universe makes BI model a prime candi- 
date for studying the possible effects of an anisotropy in the early universe on modern-day data 
observations. In view of what has been mentioned above we choose the gravitational part of the 
Lagrangian (12.21) in the form 

^-5. ("J 
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where R is the scalar curvature, K = SnG being the Einstein's gravitational constant. The gravita- 
tional field in our case is given by a Bianchi type I (BI) metric 

d s 2 = dt 2 - a 2 dx 2 - b 2 dy 2 - c 2 dz 2 , (2.4) 

with a, b, c being the functions of time t only. Here the speed of light is taken to be unity. 



B. Viscous fluid 



The influence of the viscous fluid in the evolution of the Universe is performed by means of its 
energy momentum tensor, which acts as the source of the corresponding gravitational field. The 
reason for writing Jz? v f in (12.21) is to underline that we are dealing with a self-consistent system. 
The energy momentum tensor of a viscous field has the form 

T /T(m) = ( £ + P) u ^ uV ~ P + r lS vli >m;j3 + u p-.fi - Ufiu a up. a - upu a u^ a ], (2.5) 

where 

Here £ is the energy density, p - pressure, r\ and ^ are the coefficients of shear and bulk viscosity, 
respectively. Note that the bulk and shear viscosities, r\ and ^, are both positively definite, i.e., 

T]>0, ^>0. (2.7) 

They may be either constant or function of time or energy, such as : 

r\ = \A\e a , ^H 5 !^- (2-8) 

The pressure p is connected to the energy density by means of a equation of state. In this report 
we consider the one describing a perfect fluid : 

p = Ce, Ce(0,l]- (2.9) 

Note that here £ 7^ 0, since for dust pressure, hence temperature is zero, that results in vanishing 
viscosity. 

In a comoving system of reference such that = ( 1 , 0, 0, 0) we have 

*0(m) = £ ' ( 2 - 10a ) 



I? W = V + (2.10b) 



T 3 \ m) = -p' + 2r} C -. (2.10d) 
Let us introduce the dynamical scalars such as the expansion and the shear scalar as usual 

= 4> o 2 = l -o llv o^\ (2.11) 
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where 



<V = \{un;aPv +«v;a^f) - (2.12) 
Here P is the projection operator obeying 



P 2 =P (2.13) 



For the space-time with signature (+, — , — , — ) it has the form 

PllV = g^lV —UpUv, Py=8y-U^U V . (2.14) 

For the BI metric the dynamical scalar has the form 



and 



d b c x 

= - + - + - = -, (2.15) 
a b c x 

a L b L c z 3 



C. Field equations and their solutions 

Variation of d2.ll ) with respect to metric tensor g^ v gives the Einstein's field equation. In 
account of the A-term for the BI space-time (12.41) this system of equations can be rewritten as 



! + - + !- = KZ?-A, (2.17a) 
b c be 



c a c a 



+ - + -- = kT^-A, (2.17b) 
c a c a 



a b db 



'-+u + '-u = kT 3~ A i (2-lVc) 
a b ab 

db be c d n 

+ + — = K7 o " A > ( 2 - 17d ) 
ab be c a 

where over dot means differentiation with respect to t and 7y is the energy-momentum tensor of 
a viscous fluid given above (12.101) . 



1. Expressions for the metric functions 



To write the metric functions explicitly, we define a new time dependent function x{t) 

T = abc=^g, (2.18) 

which is indeed the volume scale of the BI space-time. 

Let us now solve the Einstein equations. In account of (12.101) from (12. 17ab . (I2.17bl) . and (I2.17cl) 
one finds the following expressions for the metric functions explicitly H23I1 



a{t) = AiT 1/,3 exp 
b(t) = A 2 T 1/3 exp 
c(t) = A 3 T 1 / 3 exp 



(Bi/3) 
(B 2 /3) 
(B3/3) 



e -iK sndt - 




dt 

X 


(2.19a) 


e -2Kfr)dt - 




dt 


(2.19b) 


X 




e -2KfT]dt - 




dt 

X 


(2.19c) 
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where the constants A,'s and 5,'s obey the following relations 

AiA 2 A 3 = 1, 
Bi+B 2 +B 3 = 0. 

Thus, the metric functions are found explicitly in terms of x and viscosity. 

As one sees from (12. 19ab . d2. 19bb and (I2.19cl) . for x = t n with n > 1 the exponent tends to unity 
at large t, and the anisotropic model becomes isotropic one. 



2. Singularity analysis 

Let us now investigate the existence of singularity (singular point) of the gravitational case, 
which can be done by investigating the invariant characteristics of the space-time. In general 
relativity these invariants are composed from the curvature tensor and the metric one. In a 4D 
Riemann space-time there are 14 independent invariants. Instead of analyzing all 14 invariants, 
one can confine this study only in 3, namely the scalar curvature I\ = R, I2 = Rp, v R^ v , and the 
Kretschmann scalar I3 = R a pfivR a ^ v B, EfJ.. At any regular space-time point, these three 
invariants I\ , I2, h should be finite. Let us rewrite these invariants in detail. 

For the Bianchi I metric one finds the scalar curvature 



/a b c db be c d\ 

h =R= -2 -+T + - + -T + T- + — • 

\a b c ab be caJ 



(2.20) 



Since the Ricci tensor for the BI metric is diagonal, the invariant I2 = R^ V R^ V = R^Ry is a sum 
of squares of diagonal components of Ricci tensor, i.e., 



h 



(r°) 2 + ( R \y + (Riy + (rD 



>^2 



2\2 



>3\2 



(2.21) 



Analogically, for the Kretschmann scalar in this case we have 1$ = R^ v a gR a p V , a sum of squared 
components of all nontrivial R^v, which can be written as 



h 



ni) 2 + (* 02 o 2 ) 2 + (n 3 ) 2 + (* i2 i 2 ) 2 + (* 23 23 ) 2 + (* 31 31 



a' 



b c 



c a 



Let us now express the foregoing invariants in terms of x. From Eqs. (12.191) we have 

-2k f r\dt 



(Bi/3) J 



-dt 



Cli 

iij 



(i= 1,2,3 



X+Bie- 2K ^ dt 
3x 

3xx-2x 2 - xB ie - 2K ^ dt -6K7}xBie- 2K ^ dt +Bfe~ 4K 

9X 2 



/ rjdt 



(2.22) 



(2.23a) 

(2.23b) 
(2.23c) 
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i.e., the metric functions a,b,c and their derivatives are in functional dependence with x. From 
Eqs. (12.231 ) one can easily verify that [23] 

1 1 1 

/ 1 OC OC fn OC 

1 T 2 ' 2 T 4 ' 3 T 4 ' 

Thus we see that at any space-time point, where x = the invariants h,h, and I3 become infinity, 
hence the space-time becomes singular at this point. 



D. Equations for determining x 

In the foregoing subsection we wrote the corresponding metric functions in terms of volume 
scale x. In what follows, we write the equation for x and study it in details. 

Summation of Einstein equations (12. 17ab . (12. 17bb . d2.17cl) and 3 times d2.17dl) gives 

3 3 

x- -k^x = -K(e- p)x-3Ax. (2.24) 

For the right-hand-side of (12.241) to be a function of x only, the solution to this equation is well- 
known 12711 . 

The energy-momentum conservation law, i.e., 

^ V ; v = t^ v +t v pv tp -rj v r p v = o, (2.25) 

in our case gives the following equation for £ : 



£ + -co - (£ + U ) ^ + 477 (fcr ° - A) = 0, (2.26) 

X 3 X z 



where 

CO = e + p, (2.27) 

is the thermal function. 

Defining a generalized Hubble constant H : 

x d b c 

- = - + - + - = 3H. (2.28) 
x a b c 

the Eqs. (12.241) and (12.261) in account of (12.101) can be rewritten as 

H = ^(3%H-co)-(3H 2 -K£+A), (2.29a) 

£ = 3H(3^H-co) +4r}(3H 2 -K£+A). (2.29b) 
In terms of dynamical scalars and a the system (12.291 ) takes a very simple form 

= ^(|0_ Q) )_3 CT 2 j (2.30a) 

£ = e(£0-fi>) +47] a 2 . (2.30b) 
Note that the Eqs. (12.301 ) coincide with the ones given in Jl] 
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III. QUALITATIVE ANALYSIS AND SOME SPECIAL SOLUTIONS 

In this subsection we simultaneously solve the system of equations for T, H, and e. It is 
convenient to rewrite the Eqs. (12.281) and ( 12.291 ) as a single system : 

i = 3Ht, (3.1a) 

H = ^(3%H-co)-(3H 2 -K£+A), (3.1b) 

e = 3H(3^H-Q)) +4r}(3H 2 -K£+A). (3.1c) 
In account of (I2.27I) . (I2.8I) and (12.91 ) the Eqs. (13.11) now can be rewritten as 

t = 3Hz, (3.2a) 

H = |(35e%-(l + C)e)-(3// 2 -fce+A), (3.2b) 

e = 3#(3fie%-(l + C)e) +4Ae a (3// 2 - kc+A). (3.2c) 

The system (13.11) have been extensively studied in literature either partially H, 12, 3 or in 

general [11]. In what follows, we consider the system (13.11) for some special choices of the param- 
eters. 



A. Qualitative analysis 



Following Belinski and Khalatnikov [ 1 1] let us now study the characters of the solutions of the 
dynamical system (13.11) or (I3.2I) . We first rewrite the system (I3.ll) . namely (I3.1bl) and (B.lcl) in the 
matrix form : 

H\ = (k/2 -\\( 3%H-(Q 
£ J ~ \ 3H 47] J \ 3H 2 -K£+A 

Note that unlike the system studied by Belinski and Khalatnikov the system in consideration con- 
tains a Cosmological constant A. 



(3.3) 



1. General properties of the system 

Easy to note that the solutions cannot intersect the axis £ = 0, since e| e =o = 0, as well as the 
parabola 

3// 2 -K£+A = 0, (3.4) 

as far as (13.41 ) is itself the integral curve. Thus, starting from the point (//,£) = (+°°,0), the 
solutions cannot enter into the "prohibited region" inside the parabola (I3.4I) . Whether they may 
achieve H < depends on the value of A. Note that, unlike the system considered by Belinski et. 
al 111 ill the system in this report contains a nonzero A term. 



2. Critical points of the dynamical system 

a) By virtu of linear independence of the columns of the matrix of the Eq. (13.31) the critical 
points are the solutions of the equations 

3%H-G) = 0, (3.5a) 

3H 2 -K£+A = 0. (3.5b) 
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i.e., they necessarily lie on the parabola (I3.4I) . Solutions to the system (13.51) will be the roots of the 
equation 

3kB 2 £ 1+2 P - (1 + Q 2 e 2 - 3Afi 2 £ 2/3 = 0, (3.6a) 

H = (3.6b) 

The quantity of the positive roots of the Eq. (13.61) according to Cartesian law is equal to the 
number of changes of sign of coefficients of equations or less than that by an even number. So, for 

A < and 1/2 < j8 < 1 (Fig® Fig® 
or A>0 and j8 < 1/2 (Fig&Fig® 

the number of roots is either 2 or zero. For the remaining cases 

A<0 and /3 > 1 (Fig®, 
A<0 and /3 < 1/2 (Fig®, 
A>0 and j8 > 1/2 (Fig® 

there exists only one root. The corresponding pictures of the phase curves are given in figures 
cited above. The critical points are denoted by small circles. Note that here we consider the case 
with rj = 0, i.e., A = 0. In case if r\ ^ 0, with the increase of A the separatrix of the saddle tilts 
(inclines) to the left. Since the overall picture for A ^0 remains qualitatively unaltered, we only 
show the corresponding phase portrait for two cases, namely Fig. [8] corresponds to Fig. [2 Fig. 
|9] corresponds to Fig. |4j Note that for numerical calculations we set fc = 1, £ = 0.333 (if not 
mentioned otherwise). In the Figs. E-Qr/ is taken to be zero. Note that in the Figs. E and T stand 
for e and T, respectively. 

Since, the equation for £ only contains rj, the energy density for nontrivial 77 undergoes essen- 
tial changes, whereas H and x remain virtually unchanged. 

The types of critical points lying on the integral curve alternate: . . . saddle, attracting knot, 

saddle So it is sufficient to consider the case with maximum number of roots. Taking into 

account the Eqs. (13. lcb and (13.41) let us now calculate 

£ ,. 3B£pyjK£-A-£(l + tl) 

lim = lim 

. -(1 + C) < 0, j8<1/2, 

3BVKEW- 1 ) - Ae 2 -(1 + 0= { (3.7) 

+00 > 0, j8 > 1/2. 

So, the latest critical point for /3 < 1/2 is attracting knot and for /3 > 1/2 is saddle. 

b) It is obvious that if Lambda < the points of intersection of the boundary are the critical 
points 



H = ±V-A/3, (3.8a) 
e = 0. (3.8b) 

c) For H < there may exist critical points , if the columns of the matrix of (13.31) are linearly 
dependent. In that case the critical points are the roots of the equation 

3k(C, - l)e + 6K 2 AB£ a+ P +Sk 2 A 2 £ 2cc + 6A = 0, (3.9) 



Bianchi type I universe with viscous fluid 



9 



and 



2 

H = — kAb". 
3 



(3.10) 



In case of r\ = the roots of the characteristic equation 




(3.11) 



are 



Ml,2 



3jc£±V9k 2 £ 2 -48A(1 + 



(3.12) 



4 



The critical point (//, e) = (0, 2A/[k(1 - £)]) is of type divergent focus if A > 9?c 2 ^ 2 /[48(l + 
Q] or divergent knot if A < 9fc 2 <^ 2 /[48(l + Q] . 

In the cases illustrated in Figs. |3]and|7J H — > °° and £ — > °° as f — >■ °°, whereas, for the cases 
given in Fig. |6]one observes increasing oscillation bounded by the attracting parabola (I3.4I) . 



For A < the solutions starting from the upper half -plane H > cannot enter into the lower 
one. For A > some of the solutions may enter into the lower half-plane through the segment 
H = and < £ < A and never returns back, since H\h=o < 0. 



In this subsection solutions to the system of equations d3.ll ) has been obtained numerically. 
Evolution of the Hubble constant H, energy density £ and volume scale T corresponding to the 
cases studied above with different B, /3 and A has been illustrated in the Figs. El -El As one 
sees, for a negative A the volume scale x expands exponentially, whereas, for a positive A there 
exist solutions where % initially expands and after reaching some maximum begins to contract and 
finally collapses into a point, thus giving rise to space-time singularity. Beside this, as one sees 
from Fig. [TH a suitable choice of initial conditions gives rise to a singularity-free oscillatory mode 
of expansion of the Universe. 



Let us first consider the case when the real fluid possesses the bulk viscosity only. The corre- 
sponding system of Eqs. can then be obtained by setting r\ = in (13.11) or A = in (13 .2b . In this 
case the Eqs. (13. 1 ab and (G.lbl) remain unaltered, while (13 . 1 cb takes the form 



3. Integral curves 



B. Numerical solutions 



C. exact solutions 



In this subsection we consider some special cases allowing exact solutions. 



1. Case with bulk viscosity 



£ = 3H(3^H-0)). 



(3.13) 



In view of (13.131) the system (13.11) admits the following first integral 

T 2 ( K£ - 3H 2 - A) = Ci , Ci = const. 



(3.14) 
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The relation (13.141) can be interpreted as follows. At the initial stage of evolution the volume scale 
T tends to zero, while, the energy density £ tends to infinity. Since the Hubble constant and the A 
term are finite, the relation (13.141) is in correspondence with the current line of thinking. Let us see 
what happens as the Universe expands. It is well known that with the expansion of the Universe, 
i.e., with the increase of T, the energy density £ decreases. Suppose at some stage of expansion 
T — > oo, hence £ — > 0. Then from (13.141) follows that at the stage in question 

3H 2 + A^0. (3.15) 

In case of A = 0, we find H = 0, i.e., in absence of a A term, once T — > °°, the process of evolution 
is terminated. As one sees from (13.151) . for the H to make any sense, the A term should be negative. 
In presence of a negative A term the evolution process of the Universe never comes to a halt, it 
either expands further or begin to contract depending on the sign of H = ±\/— A/3, A < 0. 
Let us now consider the case when the bulk viscosity is inverse proportional to expansion, i.e., 

<^e=C 2 , C 2 = const. (3.16) 

Now keeping into mind that 6 = t/t = 3H, also the relations (13 .lab . (12.271) and (12.91) the Eq. 
(13.131) can be written as 

£ i 

(3.17) 



C 2 -(1 + C)e t 
From the Eq. (13.171) one finds 

£ = T -^[C 2 + C 3 T-( 1+ £)] ) (3.18) 

with C3 being some arbitrary constant. Further, inserting £ from (13.181) into (12.241) one finds the 
expression for % explicitly. 

Taking into account the equation of state (12.91) in view of (13.161) and (13.181) . the Eq. (12.241) 
admits the following solution in quadrature : 

dx 

t + t , (3.19) 



ci+cy+c\T^ 

where C 2 and to are some constants. Further we set to = 0. Here, Cq = 3kC 2 /(1 + Q — 3 A and 
C\ = 3kC3/(1 + Q. As one sees, Cq is negative for 

A>jcC 2 /(1 + 0- (3-20) 

It means that for a positive A obeying (13.201) (we assume that the constant C 2 is a positive quantity) 
T should be bound from above as well. It should be noted that for a suitable choice of C| and To 
(the initial value of t), it is possible to obtain oscillatory mode of expansion with T being always 
positive, i.e., a singularity free evolution of the Universe. The phase portrait of the (//, e) plane 
and the evolution of the BI Universe corresponding to this portrait allowing oscillatory solutions 
are given in Figs. [TOlandfTTT 

As a second example we consider the case, when £ = 1 . From (13.191) one then finds 



t(?) = (exp( V / C °0-C 2 2 exp(JC °0)/(2VO, C° > 0, (3.21a) 



T(t) = (C 2 2 /J\C° \)sm(J\C Q\t). C°q<0. (3.21b) 
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Taking into account that Cq > for any non-positive A, from (13 .21 all one sees that, in case of 
A < the Universe may be infinitely large (there is no upper bound), which is in line with the 
conclusion made above. On the other hand, Cq may be negative only for some positive value of 
A. Thus we see that a positive A can generate a oscillatory mode of expansion of a BI Universe. 
The oscillation takes place around the critical point (//, e) = (0, (2A— kC?)/\k (1 — £)]) having 
the type of cycle under the condition A > + £). It was shown in [20, 21] that in case of a 

perfect fluid a positive A always invokes oscillations in the model, whereas, in the present model 
with viscous fluid, it is the case only when A obeys (13.201) . Unlike the case with radiation where 
BI admits a singularity-free oscillatory mode of evolution, here, in case of a stiff matter one finds 
the BI Universe first expands, reaches its maximum and then contracts into a point, thus giving 
rise to space-time singularity. 

2. Case with shear and bulk viscosity 

Let us now consider the general case with the shear viscosity r\ being proportional to the ex- 
pansion, i.e., 

rjoc9 = 3H. (3.22) 

We will consider the case when 

n = ~H. (3.23) 

In this case from (B.lbl) and d3.1cl) one easily find 

3H 2 = K£ + C 4 , C 4 = const. (3.24) 

From (13.241 ) it follows that at the initial state of expansion, when e is large, the Hubble constant 
is also large and with the expansion of the Universe H decreases as does £. Inserting the relation 
d3~2"4"l) into the Eqs. (l3~Tbl one finds 



/ 



(3.25) 



Vah 2 +bh+c 

where, A = -1.5(1 + Q, B = 1.5k£, and C = 0.5C 4 (£ - 1) - A. For | being a constant, (l3~25b 
admits sinusoidal solution, i.e., H evolves oscillatory. Further, from (I3.1a|) one finds the expression 
for T, which is exponential one accompanied by a sinusoidal mode B23I1 . 



IV. CONCLUSION 

We investigated the cosmological solutions to the equations of General Relativity for the ho- 
mogeneous anisotropic Bianchi type I model by taking into account dissipative processes due 
to viscosity and Cosmological constant (A term). A detailed analysis showed that the viscosity, 
as well as the A term exhibit essential influence on the character of the solutions. The classifi- 
cation of the solutions was pursued for the viscosity being some power law of energy density, 
namely, rj =Ae a and E, = Be^ , It was noticed that for A < the Universe expands forever with 
a logarithmic velocity H, which, depending on the viscosity either becomes constant or increases 
infinitely. In the process behavior of the energy density e is analogous to that of H except the case 
when e — > 0. For A > 0, beside the variants mentioned above, there exists few other possibilities: 
contraction of the Universe into a point, thus giving rise to a space-time singularity; a regime of 
increasing oscillation corresponding to suitable initial conditions. It was also noticed that a special 
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case with A > 0, rj = and £,H = const, the model admits a singularity-free oscillatory mode of 
expansion. 
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FIG. 1 : Phase diagram on H — e plane for 
jS = 1.5, A = -.933, fi = .720. 
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FIG. 2: Phase diagram on H — e plane for 
jS = .75, A = -.707, £ = .589. 
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FIG. 3: Phase diagram on H — e plane for 
J3 = .75, A = -.707, fi = .667 



FIG. 4: Phase diagram on H — e plane for 
j3 = .05, A = -.785, B = .451. 




FIG. 5: Phase diagram on H — e plane for FIG. 6: Phase diagram on H — e plane for 

j8 = .05, A = .317, B = 0.933 J8 = .05, A = .317, B = .667. 
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FIG. 9: Phase diagram on H — e plane for 
j3 = .05, A = -.785, B = .451,A = 1, cc = 

1, K=l. 
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FIG. 10: Phase diagram on H — e plane for 
A = 3, £ = 0.333, C 2 = 1, C 3 = 1. 
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FIG. 11: Evolution of the BI Universe corresponding to the phase diagram 
given in Fig. El| As one sees, the BI Universe in this case undergoes an 
oscillatory mode of expansion 
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FIG. 12: Evolution of the FIG. 13: Evolution of the en- FIG. 14: Evolution of the 

Hubble constant H with pa- ergy density e with parame- volume scale X with param- 

rameters as in Fig. [I] ters as in Fig. [I] eters as in Fig. [I] 




FIG. 15: Evolution of the FIG. 16: Evolution of the en- FIG. 17: Evolution of the 

Hubble constant H with pa- ergy density e with parame- volume scale X with param- 

rameters as in Fig. El ters as in Fig. eters as in Fig. |2] 




FIG. 18: Evolution of the FIG. 19: Evolution of the en- FIG. 20: Evolution of the 

Hubble constant H with pa- ergy density e with parame- volume scale x with param- 

rameters as in Fig. |3] ters as in Fig. |3] eters as in Fig. |3] 




FIG. 21: Evolution of the FIG. 22: Evolution of the en- FIG. 23: Evolution of the 

Hubble constant H with pa- ergy density e with parame- volume scale X with param- 

rameters as in Fig. |4] ters as in Fig. 0] eters as in Fig. |4] 
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FIG. 24: Evolution of the 
Hubble constant H with pa- 
rameters as in Fig. [5] 




FIG. 25: Evolution of the en- 
ergy density e with parame- 
ters as in Fig. |5] 
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FIG. 26: Evolution of the 
volume scale X with param- 
eters as in Fig. |5] 
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FIG. 27: Evolution of the 
Hubble constant H with pa- 
rameters as in Fig. |6] 




FIG. 28: Evolution of the en- 
ergy density e with parame- 
ters as in Fig. |6] 




FIG. 29: Evolution of the 
volume scale % with param- 
eters as in Fig. |6] 




FIG. 30: Evolution of the 
Hubble constant H with pa- 
rameters as in Fig. 




FIG. 3 1 : Evolution of the en- 
ergy density e with parame- 
ters as in Fig. [7] 




FIG. 32: Evolution of the 
volume scale t with param- 
eters as in Fig. 



